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Abstract 



Unlike isolated /3-cells, which usually produce continuous spikes or fast and 
irregular bursts, electrically coupled /3-cells are apt to exhibit robust burst- 
ing action potentials. We consider the noise induced by thermal fluctuations 
as well as that by channel gating stochasticity and examine its effects on the 
action potential behavior of the /3-cell model. It is observed numerically 
that such noise in general helps single cells to produce a variety of electrical 
activities. In addition, we also probe coupling via gap junctions between 
neighboring cells, with heterogeneity induced by noise, to find that it en- 
hances regular bursts. 

Key words: Thermal fluctuation; Channel gating stochasticity; Hetero- 
geneity; Gap junction 
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Introduction 



Bursting action potentials, which are characterized by rapid firing inter- 
spersed with quiescent periods in pancreatic /3-cells, play a central role in 
the secretion of insulin, the hormone for glucose homeostasis. It has been 
reported that isolated /3-cells actually show continuous spikes or fast and ir- 
regular bursts while /3-cells in a cluster or in an intact islet produce 
regular bursting action potentials 

Mm 

. As for the correlations be- 
tween the electrical activity on the cell membrane and insulin secretion (jj), 
the robust bursts appear more effective in maintaining glucose homeosta- 



sis than continuous spikes, since cou 



better than isolated /3-cells 



(9, 



lid 



jled /3-cells can control insulin release 



Ul). 



However, the question as to whether bursting is an endogenous prop- 
erty of individual /3-cells or of a cluster still remains to be answered, which 
has attracted a number of investigations. Among proposed explanations is 
the channel-sharing hypothesis, which postulates that current fluctuations 
arising from channel gating stochasticity prevent single cells, originally ca- 
pable of bursting, from bursting, but when they are electrically coupled, 
the perturbin g eff ects are shared by neighbors and the regular bursting is 



recovered 



(3, 3,0). 



noise, recent research 



In contrast to this hypothesis of negative effects of 
12, Q, 3) has established that noise can 



£3, Q, 



>lay a constructive role in many biological systems including /3-cell bursting 
2w . The heterogeneity hypothesis, providing another explanation, was also 
postulated by the same group. According to it, when heterogeneous cells, 
each of which produces continuous spikes or bursts depending upon such cell 
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parameters as the size, channel density, etc., are coupled, those cells in the 
cluster exhibit more pronounced bursts. This gives a useful insight into the 
functioning of heterogeneous cell populations (21). 

In this study, we expand the concept of heterogeneity and probe how 
such general heterogeneity enhances bursting. It is proposed that noise in- 
duces heterogeneity in otherwise homogeneous individual /3-cells, which in 
turn assists the /3-cells to produce robust bursts when they are coupled. 
Existing studies have mostly focused on the synchronizing role of coupling 
22I 23L the slow dynamics, which has a period about 10 to 60 seconds, is 
synchronized successfully between adjacent cells. In contrast, we focus here 
on the fact that rapid firing in the active phase of bursting is asynchronous 
between neighbors (24) and these fluctuating currents through the gap junc- 
tion act like noise, enhancing the robust bursting action potential. It is also 
presented that various action potentials of single /3-cells are embodied with 
optimal noise induced by thermal fluctuations or by ionic channel gating 
stochasticity. In particular, noise stimulates occasionally itself to produce 
fast bursts in a single cell. 

There are four sections in this paper: In the second section the math- 
ematical model for /3-cells is introduced and the simulation method is de- 
scribed. The third section is devoted to the effects of random noise in cur- 
rents and of voltage-dependent noise in single cells while the fourth section 
examines how coupling between cells influences the electrical activity of a 
cell. Finally, main results are summarized and discussed in the last section. 
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Model and Methods 



Mathematical model for a /3-cell 



As the Hodgkin-Huxley model (|25T ) describes the electrical activity on the 



cell membrane with ion channels, a few mathematical models for /3-cells, 

nnQ 

( 2n . 1271 . 1281 ) of the ion channels in 



based on the electrophysiological data 



/3-cells, have been 



proposed. Although there are simple models using two- 



dimensional maps (|2£ , 



30 



311 ) , we consider the Sherman model, which allows 



direct physical interpretation 



20 



32) 



The model is described by the current balance equation between capac- 
itive and ionic currents: 



dV 



-ICa(V) ~ I K (V, N) - I KiA TP)(V, P) 

-Is(V,S), 



where Cm and V denote the membrane capacitance and the membrane 
potential, respectively. The activation variable N and the slow variable S 
are governed by 

dN 

rjv— = N^-N 
dt 

rs^ = S^-S (2) 

with appropriate relaxation times tm and rg, which are taken to be constants 
for simplicity. The fraction P of open K(ATP) channels may also be regarded 
as a constant for the moment [see Eq. Ionic currents here are fast 
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voltage-dependent L-type Ca + current Ic a , delayed-rectifier K + current 
Ik, ATP-blockable K + current Ik(ATP), and very slow inhibitory potassium 
current Is' 



I C a(V) = gCaM^VW - Vca) 

I K (V,N) = g K N(V-V K ) 
Ik(atp)(V,P) = 9k(atp)P{V ~V K ) (3) 
Is(V,S) = g s S(V-V K ). 

Ica and Ik are responsible for generating action potentials; Ic a is assumed 
to respond instantaneously to a change in the membrane potential, whereas 
Ik is governed by the dynamics of the activation variable N via Eq. |2J 
Ir(ATP) is the background current with voltage-independent conductance 
9k (ATP)] this determines the plateau fraction, i.e., the ratio of the active 
phase duration to the burst period. For example, as 9k(ATP) decreases 
under high glucose concentration, there are only active phases without silent 
phases. Is is a phenomenological current representing slow dynamics in 
the bursting action potential. This model thus assumes that single /3-cells 
originally contain the slow dynamics, which works just under the appropriate 
condition. Biological candidates for such slow dynamics include slow free 

Q). 



Ca + dynamics 



33j) and ATP metabolism (|34f). Finally, M^, Noo, and S 



of the voltage-dependent activation are defined to be 



Xoo{v) ~ i + ex P [(v x -v)/e x y (4) 
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where X denotes M, N, or S. 

This set of coupled nonlinear differential equations in Eqs. H - El has 
been analyzed in detail There it is noted that S responds on a 

much slower time scale than V and N because ts has the time scale of 
several seconds compared with the milli-second time scale in firing. Then 
S is regarded just as a parameter, and the dynamics of the fast subsystem 
on the two-dimensional phase space of V and N is analyzed. Furthermore, 
after eliminating one degree of freedom by substituting iVoo to N, the whole 
behavior of this model may be analyzed approximately with fast variable V 
and slow variable S. 

Numerical details 

Integration of differential equations including noise demands some caution, 
and is commonly achieved via the Euler method. For better efficiency, we 
employ the Euler method for integrating the noise term, combined with the 
second-order Runge-Kutta method for other terms. In order to be concrete, 
we consider the one- variable problem 



fin- 

- = /(*)+£(*), (5) 



where f(x(t)) is a (nonlinear) function of x, the variable of concern, and 
is the white noise with zero mean and delta- function correlations 



«(*)) = o, 

mat')) = 2D6(t-t'). 



(6) 
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Taking the time step of size At, we obtain from the equation of motion the 
value of x at time t + At: 



x(t + At) = x(t) + f ( x ®\ + /(x) At + £(t)At, (7) 



(Q). 



where x = x (t) + /(x(t))At + £(t)At (37). Although there is no gurantee 
that this algorithm should converge in general,it works fine here since the 
noise term does not depend on the variable x (|38l ). 

The white noise £ of variance D is produced by the gaussian random 
numbers with the variance a 2 determined by 



where the Dirac delta function has been represented by At 1 within the 



numerical accuracy. We thus have the relation a = w2D/At. 

In our simulations, we take At = 1 ms, which turns out to be small 
enough, and integrate the set of equations for current balance. This gives 
the time evolution of the action potential, from which the power spectrum 
is computed through the use of the fast Fourier transform technique. 



Results and Discussion 

Noise effects 

Before explaining the coupling effects, we first probe the role of noise, either 
the usual (additive) random noise or the (multiplicative) voltage-dependent 
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one. Comparison of the effects of such noise helps us to understand better 
the coupling effects. 

Random noise 

Among many kinds of noise on the cell membrane, the simplest case is the 
random noise, which may come from thermal fluctuations (see below). When 
such random current fluctuations are present on the membrane, the current 
balance equation in Eq. ^ is generalized to 

Cm 1E~ = -W^,iV,.S)-£(t), (9) 

where Ii on represents all the ionic currents on the right-hand side in Eq. ^ 
and the noise current satisfies Eq. |f)]with the variance denoted by D^. 

Figure Q exhibits the solution of the set of coupled differential equations 
in Eqs. [2] and under various strengths of the random noise. It is observed 
that single /3-cells produce various electrical activities according to the value 
of ttv in Eq. [21 which lies in the narrow range 4 to 11 ms depending on the 
membrane potential (27). When the time constant tn of delayed-rectifier K + 
channel activity exceeds 11.0 ms, the /3-cell produces regular spiking action 
potentials in Fig. [I\A, while for tm below 10.0 ms faster repolarization does 
not allow enough time for the slow variable S to decrease, yielding bursting 
action potentials [see Fig. ^ C]. In the intermediate regime of T7v=10.2ms, 
Fig. ^ B shows that spiking action potentials are generated but the bursting 
property is resident. As an appropriate amount of noise comes into play, in 
particular, the regular spikes in Figs. ^ A and B and bursts in C change 
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into fast bursts in E, irregular spikes in G, or irregular bursts in H and /. 

To explain these phenomena, we note two thresholds of the slow variable 
S: One is the upper threshold above which the membrane potential is falling 
into the resting potential; the other is the lower threshold above which the 
membrane potential begins to fire. At the moment that fluctuations take 
negative values, they may assist the repolarizing membrane potential to re- 
main above the lower threshold before the membrane repolarizes completely 
and then, depolarizes slowly to the lower threshold. This induces occasion- 
ally consecutive firing in Fig. Q G or even fast bursts in Fig. ^ E for the 
/3-cell in the critical parameter range, i.e., T7v=10.2ms. Such consecutive 
firing raises the average membrane potential for a while, compared with the 
case of regular spikes. Hence the value of Soo becomes large, and conse- 
quently S grows with the delay represented by the time constant rg. When 
it goes over the upper threshold, the membrane potential returns to the rest- 
ing potential. At the same time, now becomes small and S reduces to the 
lower threshold. During this period of S varying from the upper threshold 
to the lower one, the membrane potential stays in the silent phase. When S 
comes to the lower threshold, the membrane potential starts to depolarize 
and fire. Repetition of these processes simply constitutes the fast bursts. 
As the noise level is raised further, the slow variable S may start to increase 
before it reaches the lower threshold, assisted by the fluctuations taking neg- 
ative values. Similarly it may start to decrease before it reaches the upper 
threshold due to positive fluctuations. In consequence, irregular bursts in 
Figs-E-^ an d / can thus be induced. When fluctuations become sufficiently 
strong and dominant, such a role of noise, turning on the slow dynamics of 
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S, is concealed and the membrane potential appears noisy. Here it is notable 
that under optimal fluctuations, there exists the critical parameter range in 
which the difference between the upper and lower thresholds is small and 
the dramatic effect of fast bursts is produced; similar results were obtained 
in a recent study (Q). 

It is revealing to examine the power spectra of the obtained action po- 
tentials, computed through the use of the fast Fourier transform technique 
for various noise levels and displayed in Fig. |2J In particular, Figure |2] B 
manifests that the regular spiking action potential of frequency 2 Hz in the 
absence of noise has changed into fast bursts containing oscillations of 0.2 Hz 
and 5 Hz at moderate noise levels. 

To characterize the positive/negative role of noise in bursting, we define 
the bursting tendency according to B = log[P(f b)/T , (0)], where V(/b) is 
the power spectrum at the bursting frequency fs and V(0) is the background 
intensity at Hz. Figure 01 shows the behavior of the bursting tendency B 
with the noise level, manifesting the noise effects on bursting. 

Finally, one may ask whether thermal fluctuations known to generate 
white noise are enough to induce the fast bursts, irregular bursts or spikes, 
observed in our simulations. In simulations, the variance is taken in 
the range 10 -29 J/O ~ 10~ 27 J/O. In reality, noise currents due to thermal 
fluctuations can be estimated via the fluctuation-dissipation theorem: = 
ksT / R. This gives ~ 10 -29 J/O when R is taken to be a few giga ohms 
(GO) or less. Accordingly, thermal fluctuations alone may not be enough 
to induce irregular spikes or bursts. Nevertheless, it appears possible that 
thermal fluctuations actually expedite the emergence of fast bursts when 
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the cell lies in the critical parameter regime. 
Voltage-dependent noise 

As another simple type of noise, one can consider the voltage-dependent 
fluctuations, which are closely related to the channel gating stochasticity 
(see below). In the presence of such multiplicative noise, the current balance 
condition in Eq. Stakes the form 



where Ii on also represents all the ionic currents in Eq. ^ and n(t) is the 
Gaussian white noise, again satisfying Eq. with variance D v . Solving 
numerically the coupled differential equations given by Eqs. El and at 
various noise levels with r^r set equal to 11 ms, we obtain the results, which 
are illustrated in Fig. HJ Note the overall similarity to the case of random 
(additive) noise shown in Figs.^-D and G. 

When the voltage-dependent noise stimulates the cell membrane, irreg- 
ular spikes arise, similarly to the case of random noise, if its amplitude 
multiplied by the voltage difference (V— Vk) is comparable to the ampli- 
tude of random noise. In fact, voltage-dependent noise may be regarded 
simply as the noise weighted more in the active phase of the membrane 
potential than in the silent phase. When taking negative values, there- 
fore, fluctuations boost firing more effectively in the active phase and con- 
tribute less to the erratic evolution of the resting potential in the silent 
phase. Such voltage-dependent (multiplicative) noise may arise from ion 




I ion (V,N,S)- V (t)(V -V K ) 



(10) 
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channel gating stochasticity, since currents through channels depend upon 
the membrane potential difference. If the number of channels is sufficiently 
larg e, the channel stochasticity can be described by a Langevin equation 
( 3pJ . |4(J . l4lf ). Specifically, the stochasticity of K(ATP) channels has been 
considered (j^j). In the expression for the ATP-dependent K + current, 
Ik (ATP) = 9K(ATP)I 3 (y ~ Vk), the opening ratio P, which is no more con- 
stant, evolves according to 



21(1 _P)_^ P + (11) 
at Tp Tp 



where 71 /rp and 72/ip represent the rates for a closed channel to switch 
to the open state and vice versa, respectively. Note that 71 and 72 thus 
determine the equilibrium ratio between the open state and the closed one. 
Fluctuations in the opening ratio are described by the Gaussian white noise 
£(i) satisfying Eq. H^with the variance 



D _ 71(1 - P) +72P 7172 ^ 

2 TP N K(ATP) T P N K(ATP){ll+l2Y 



where Nk(atp) is the total number of ATP-dependent K + channels in a 
/3-cell Id). 

Solving Eq. I1H we obtain that P fluctuates around the equilibrium value 
Po, taken to be 0.5 in our simulations: P(t) = Po + fj(t). Here fj(t) is colored 
noise, characterized by the variance 

WW)) = D^e-^ - je-^ +t ')} (13) 
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with 7 = (71 + 72)/t p and = D^/j 2 (see Appendix for details). Note 
that the firing time scale is comparable to the correlation time 7 -1 of the 
noise fj(t) (see Fig. EJ. Consequently this colored noise is more effective to 
induce several consecutive firings, which resemble irregular burst, than the 
white noise. In particular, the modules of several spikes are observed to 
become longer as the correlation time 7 -1 is increased. Figure |H1 shows the 
behaviors in the presence of the channel-gating noise £ (t) for two different 
channel numbers. In this case of multiplicative colored noise, modules of 
spikes arise more efficiently than in the case of mutiplicative white noise 
shown in Fig. |1J Further, it is also found that stronger gating fluctuations 
from less channels (iVjr(ATP) = 500) in Fig. El B give rise to modules of more 
rapid spikes, compared with the case N K ^ ATP ^=2500 in Fig. 0^4. 

Similar results can be obtained with fluctuations in the Ca 2+ channels 
and in the delayed-rectifier K + channels although they act somewhat dif- 
ferently from the fluctuations in the ATP-blockable K + channels (data not 
shown) . 

It is thus concluded that noise generates diverse firing patterns in single 
/3-cells. In a real (physiological) islet, however, /3-cells are not isolated but 
coupled with each other, making it desirable to consider coupled /3-cells and 
to investigate effects of noise together with those of coupling. This will be 
the subject of the next section. 
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Coupling effects 

We consider two cells coupled with each other via a gap junction. With the 
coupling incorporated, Eq. ^is extended to the coupled equations: 



Cm ^T = - I ion(Vi,N 1 ,S 1 ,P 1 )-g c (V 1 -V 2 ) 

CM ^t~ = - I ion(V2,N 2 ,S 2 ,P2)-9c(V2-V 1 ), (14) 



where the subscripts 1 and 2 are the cell indices, ij on again denotes all the 
ionic currents, and gc is the coupling conductance. Note that the hetero- 
geneity between both cells is accommodated in the K(ATP) channel opening 
ratio P. Namely, the noise associated with channel gating stochasticity in- 
duces continuously heterogeneity between the cells. 

We thus have eight coupled differential equations, which consist of Eqs.|2] 
and ^2 for each cell and Eq. 1141 for eight variables (V, N, S, and P for each 
cell). Integration of these coupled equations yields the results displayed in 
Fig. for the channel-gating noise of variance = 4 x 10~ 4 s _1 given by 
Eq.^Jand for three values of the coupling conductance: gc = 50 pS, 110 pS, 
and 200 pS. Revealed is the optimal coupling strength for longer bursting 
periods: While weak coupling is not enough to couple individual cells and 
to generate consecutive firin g, to o strong coupling tends to make the cluster 
behave as a single large cell (22). 

Robust bursts emerge as a consequence of the competition between het- 
erogeneity and coupling J23L On one hand, the coupling term in Ea, ll4l helps 
the two cells to act synchronously; on the other hand, it also plays the role 
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of stimulating noise, which acts strongly on the two cells with asynchronous 
phases. The perfect asynchrony results from the harmony of coupling to be 
similar and heterogeneity to be different (see Fig. IB}. Namely, the coupling 
currents between asynchronous neighboring cells give rise to consecutive fir- 
ing; this in turn increases the upper threshold of the slow variable S above 
which firing disappears. As S grows up toward the increased upper thresh- 
old, it takes longer to reduce down to the lower threshold. This larger rising 
and falling divides more clearly the active and silent phases in the membrane 
potential, and accordingly induces robust bursting action potentials with pe- 
riods longer than 20 s. Note that in the absence of coupling we have not been 
able to observe bursting periods longer than 10 s (see Figs. HUH) (Parameter 
values different from those in Table I may yield bursting periods somewhat 
longer than 10 s even in a single cell. In this case, the coupling gives rise 
to robust bursting of even longer periods, say, 30 s, still demonstrating its 
crucial role in generating regular bursts.) 

In the two-cell model here the optimal value of the coupling conduc- 
tance is observed to be gc = 110 pS. As the number of cells is increased, 
however, more heterogeneity is introduced, which should be matched by 
stronger coupling to generate robust bursts with longer periods. Although 
the detailed investigation is beyond our computing capacity, we have per- 
formed multi-cell simulations, which indeed confirms such an increase of the 
optimal coupling conductance. For example, the optimal conductance in the 
system of 1000 cells turns out to be 100 to 300 pS (data not shown), which 
coincides with experimental results of the gap junctional conductance j^ ). 

These features of the coupled cells do not change much in the presence 
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of the voltage-dependent noise instead of the channel-gating noise, except 
that the channel-gating noise is more efficient for robust bursting than the 
voltage-dependent one, as shown in Fig.Elfor D v = 10" 24 J/ft-V 2 and D| = 
0. Note also that the coupled cells depicted in Figs. [7| and EJ] do not burst 
in the absence of noise-induced heterogeneity. 

Recall that in the emergence of robust bursts, the asynchrony from the 
heterogeneity induced by noise plavs an important role, which has also been 
addressed in a very recent study ( 43J) • Similar to such noise- induced hetero- 
geneity, the cell-to-cell heterogeneity associated with variations of the cell 
parameters among the cells is also expected to play for robust bursts (|21T ) . To 
check this, we allowed variations of the membrane capacitance Cm related 
to the cell size as well as of the channel conductance gK(ATP) an d examine 
the resulting behavior: Shown in Fig. ^1 A and B are bursts generated in 
the case of 20% variation of Cm (5.0 pF, 6.3 pF) and in the case of 10% 
variation of Qk{ATP) (1000 pS, 1100 pS), respectively. Specifically, a spiking 
cell (with Cm = 6.3 pF) is coupled with a bursting cell (with Cm = 5.0 pF) 
in Fig. El] ^4, which results in that both cells are bursting synchronously with 
a longer bursting period than that of a single cell (5.0 pF). In Fig. ^1-3, ° n 
the other hand, two spiking cells (with gK(ATP) = 1000 pS and 1100 pS) are 
coupled with each other, and both are bursting. Therefore heterogeneity is 
in general important for bursting in coupled cells, no matter whether it is 
cell-to-cell heterogeneity or induced by noise. 
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Conclusions 

We have probed whether noise and coupling serve as an appropriate stim- 
ulus for inducing the bursting action potential in pancreatic /3-cells, and 
found that they effectively call into action the inherent slow dynamics in 
individual cells. Fast bursts, irregular spikes or bursts in single /3-cells have 
been observed as the results of the noise effects. In particular the emergence 



of regular bursts assisted by an appropriate amount of noise 



(3, 



16 



17 



see Figs. [I] 



18 



ljj). In 



E and |2] B] is reminiscent of coherence resonance 
view of physiology, the consecutive firing induced by fluctuations gives rise 
to relative depolarization for a while, which is followed by the activation 
of the slow potassium channel lasting until the slow variable reaches the 
upper threshold. At this time the slow K + channel opens fully, and the 
outflux of cytosolic potassium ions gets very large, thus hindering depolar- 
ization. Accordingly, the membrane potential is compelled to stay in the 
silent phase, and the slow K + channel in turn starts to be inactivated. In 
consequence, the membrane can become depolarized as the outflux of K + 
ions reduces. Finally, firing occurs again, and consecutive firing also happens 
by the help of appropriate stimulation. As candidates for the stimulus, both 
the (additive) random noise coming from fluctuating currents and the (mul- 
tiplicative) voltage-dependent noise from the channel gating stochasticity 
have been considered. 

In particular, coupling between cells has turned out essential for attain- 
ing regular bursts with longer periods compared with the fast bursts. The 
coupling term, proportional to the potential difference between two cells, 
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operates in a similar manner to the voltage-dependent noise: It increases 
with the potential difference and thus becomes large for the cells in active 
phases, stimulating the cells like noise. On the other hand, it is small for 
perfectly synchronized cells in silent phases. The coupling also increases the 
upper threshold of S and induces robust regular bursts. 

In the analysis, the heterogeneity has been found to play an important 
role in inducing strong fluctuations during active phases, which may cause 
robust bursts. Namely, bursting in general results from the interplay of 
coupling and heterogeneity. This allows us to interpret the fact that large 



22J): As- 



cell clusters (up to the critical size) show more regular bursts (2C, 
suming a cubic islet, we have considered /3-cells arranged into an L 3 cube, 
under free boundary conditions. Adopting physiological gap junction con- 
ductance, gc = 200 pS (42), we have found that the bursting period and 
duration first increases with the size L but tends to saturate beyond L = 5 
(data not shown). Such saturation behavior may be explained as follows: 
Via the coupling through gap junctions, the number of nearest neighbors 
in the three-dimensional space is limited, e.g., to six or so; this suggests 
that the cluster above some critical size can get no more advantage of the 
heterogeneity from neighboring cells through given coupling strength. 

The Langerhans islet, however, consists of several endocrine cells in ad- 
dition to /3-cells. Other endocrine cells in an islet have been studied recently 



(44, 



451 ). and it will be of interest to study the coupling effects between origi- 



nally different a-, /?-, and <5-cells, coupled via hormones or neurotransmitters 
(3). This might give a clue to understanding the size of a Langerhans islet 
in the pancreas, which is left for further study. 
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Appendix 

Equation can be solved to give the time evolution of the opening ratio 
P: 

P(t) = P + [P(0) - PoK 7 ' + f* e-^-^mdt' 

Jo 

with Pq = 71/(71 +72) and 7 = (71 +72)/^, where -P(O) is the initial value 
of P. After sufficiently long time, we thus have P fluctuating around the 
equilibrium ratio Pq: P(t) = Pq + fj(t), where the noise fj(t) is given by 

fj(t) = f e-^-^i^dt'. 
Jo 

From the above definition of the noise fj(t), it is straightforward to derive 
its characteristics: 

ft rt' 

(#)#')> = / dre^-V / dT'e^'-^^iT)^')) 
Jo Jo 

= 2D*T7<* + *') fdre 2 ^, 
Jo 

where we have used the relation (£(t)£(t')) = 2D^5(t — t') and i denotes 
the smaller one between t and t'. We thus obtain the correlations of the 
noise f) at different times 

(fj(t)v(t')) = A 7 [7e~ 7|t ~ t ' 1 -7e- 7(t+t,) ] 



with Dfj = Z^/7 2 , which manifests the colored nature. 
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Table 1: Standard parameter values 



C M = 6.3 pF 
g K = 4000 pS 
g s = 3000 pS 
V Ca = 25mV 
V M = -20 mV 
V N = -17 mV 
V s = -22 mV 

T N = 1.1 X 10~ 2 S 

rp = 0.50 s 
7i = 1 



g Ca = 3000 pS 
9k{atp) = 1000 pS 
5c = H0pS 
1^ = -75 mV 
# M = 12 mV 
9 N = 5.6 mV 
6*5 = 8.0 mV 
t s = 20 s 
Nk(atp) = 2500 
72 = 1 
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Figure Legends 

Figure ^ 

Action potential V and slow channel activity S in single /3-cells at the noise 
level Z?£ = 0, 10 -29 , and 10~ 27 J/S7 under several values of time constant 
Tjv of delayed-rectifier K + channel activity N. All simulations have been 
performed under the standard parameter values in Table I except r^r, the 
values of which are given above. 

Figure El 

Power spectra of the action potentials for the random noise levels in Fig. ^ 
The time constant tn of the activation variable N is (A) 11.0 ms and (B) 
10.2 ms. Observed in the power spectra are main peaks together with their 
harmonics. The peak at 1 Hz, indicated by the asterisk in (B), reflects the 
tendency to form dimerization of spikes. Each power spectrum has been 
obtained from the average over 1000 samples, each having a time sequence 
of 132 seconds. 

Figure 

Bursting tendency B of /3-cells versus the noise level for several values of tat, 
corresponding to different firing patterns in the absence of noise. 

Figure |1J 

Action potential V and slow channel activity S in single /3-cells at two values 
of the voltage-dependent noise. Again parameter values in Table I have been 
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used. 
Figure 

Correlations between the action potential and multiplicative colored noise 
due to channel-gating stochasticity. The correlation time 7 -1 is taken to 
be (A) 25 ms, (B) 250 ms, and (C) 2500 ms. Note that each figure has a 
different time scale. Their corresponding power spectra are shown in (D). 
Parameter values in Table I have been used except rp. 

Figure [|0 

Action potential V and slow channel activity S in single /3-cells at two values 
of the channel gating stochasticity: (A) and (B) correspond to the channel 
number N K ^ ATP ^ = 2500 and 500, respectively. Note that Figs.El-B and|H1^4 
represent the same sample path, but with different variables plotted. Other 
parameter values have been taken from Table I. 

Figure [7J 

Figure |H1 Enlarged view of the interval between 10s to lis in Fig. Q B, 
disclosing the detailed behavior of the two membrane potentials V\ {solid 
line) and V2 {dashed line). 

Figure EH 

Bursting action potential induced by cell coupling via the gap junction of 
conductance gc = 110 pS under the voltage-dependent noise of strength 
D v = 10 -24 J/J7-V 2 . Parameter values in Table I have been used. 
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Figure I10L 

Bursting action potential induced by cell coupling, with the cell-to-cell het- 
erogeneity due to variations of the membrane capacitance Cm and of the 
ATP-blockable K + channel conductance 9k(ATP)'- (^) 20% variation of Cm 
(5.0 pF, 6.3 pF); {B) 10% variation of 9k {ATP) (1000 pS, 1100 pS). Other 
parameter values have been taken from Table I. 
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Figure 2: 
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Figure 5: 
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Figure 7: 
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